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1 Introduction 



The instantaneous Coulomb interaction in QED is well known as a fixed experimental fact. 
Instantaneous interaction is one of the indispensable attributes of the Dirac Hamiltonian ap- 
proach to QED [HE], where a variable without the canonical momentum is converted into a 
Coulomb "potential". Remember that this Hamiltonian approach was used pi E] as the foun- 
dation of all heuristic methods of quantization of gauge theories, including the Faddeev-Popov 
(FP) method [E] and all its versions used now for description of Standard Model of elementary 
particles. 

However, a contemporary reader could not find instantaneous interactions in the accepted 
version of Standard Model (SM) of elementary particles based on the FP method [6j. 

Why did the accepted version of SM omit instantaneous interactions? 

What are the fundamental principles of foundation of instantaneous interactions? 

What are physical results following from the instantaneous interactions omitted by the 
accepted version of SM? 

Responses to these questions are the topics of the present paper. 

In Section 2, we recall the Dirac Hamiltonian approach to gauge theories, where the instan- 
taneous interactions was introduced by Dirac in agreement with the priority of quantum 
principles and experimental facts. Section 3 is devoted to the statement of the problem of 
description of instantaneous interactions in SM. In Section 4, some experimental tests are dis- 
cussed as evidences of instantaneous interactions. In Appendix A, S-matrix for bound states in 
electrodynamics is derived in terms of bilocal fields. In Appendix B, the problem of the chiral 
hadronization of QCD in the spirit of the non-Abelian generalization of the Dirac approach to 
QED is discussed. 



2 Status of instantaneous interaction in gauge theories 
2.1 The frame- depended variational principle 

In order to explain the origin of the instantaneous interaction in SM, let us consider the Dirac 
approach [I] to QED given by the action 

Sqed = J d 4 x£QED, (1) 

with Lagrangian 

£qed = ~\ M„ - d u A,} 2 + $[ift - 77# - A^. (2) 

where ft = (9^7^, A^ is a vector potential, if) is the electron-positron field described by the Dirac 
bispinor, j M = eip^^ip is the charge current, and e is the electron charge. 
The action (JJJ is invariant with respect to gauge transformations [3 [E] 

A* = A, + d,\, ^ A = e +ie V, ^ A = e- JeA V>. (3) 
Physical solutions of the system of variational equations 

= dv _ qh A vi -j» = (4) 
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are obtained in a specific inertial frame of reference to initial data 1 distinguished by a unit 
timelike vector n M (ril = 1). This vector splits gauge field into the timelike A = A^n^ and 
spacelike A^ = A v — n u (A IM n fJj ) components 2 . 

The time component of Eq. ((H) with respect to A is the Gauss constraint 

AA - d d k A k = -jo, (5) 

where A = d? j Q — eip^oip is the current. In accordance with the theory of differential 
equations, the field components A cannot be a degree of freedom, because the canonical 
momentum P = 9C/d(d Ao) = is equal to zero. An exact solution of Eq. © 

a i \ 1 \a a -i- 1 f [9k A k (x ,y k ) -jo(x ,yk}] 

A (x ,x k ) = —[d k A k - jo = - — / d 6 y — . (6) 

A 4tt J \x — y\ 

is treated as the Coulomb potential field leading to the instantaneous interaction. The source of 
this potential field can be only an electric current jo- The linear terms dod k A k in this solution 
(@) cannot be considered as a physical source of the Coulomb potential. Therefore, Dirac [T] 
supposed to remove this term from the Gauss constraint JSJ) by gauge transformations 

A^ d \A) = A fl + ^A (rad) (A), ^ (rad) (^, A) = e jeA(rad)(A V, (7) 

where 

A (rad) (^) = -x* A ^^l dS dkMx °\ yk) . (8) 
A 4tt J \x — y\ 

2.2 Radiation variables as Dirac's gauge-invariant observables 

This transformation (0) determines new variables 

4 rad) (A) = A + d A {iad \A) = A -d ^d k A k , (9) 
A[ rad) (A) = A l + d l A^ d \A) = A l -d l ^d k A k , (10) 
^ (rad) (^) = e ieA(rad)(A V- (11) 

Thus, the frame-fixing A^ = (A ,A k ), the treatment of A as a classical field, and the Dirac 
diagonalization of the Gauss law [1\ lead to variables as the functionals that are invariant 
with respect to gauge transformations Q of the initial fields A^ip 

A { ™ A \A + d\) = A { ™ d) (A), ^ (rad) (^e ieA , A + d\) = ^ (rad) (^, A). (12) 



1 The physical concept of a frame of reference to initial data is defined as a three-dimensional coordinate basis 
with a watch and a ruler (and other physical devices) for measurement of time and distance (velocity, mass, 
and other physical quantities), i.e. the initial data required for unambiguous resolving equations of motion J3J|. 
Inertial means that this coordinate basis is connected with a heavy physical body moving without influences of 
any external forces. Inertial frame of reference in Minkowskian space-time is associated with the unit time-axis 

= (-^=p:, " _p ), where v is the velocity. The frame of reference — (1, 0, 0, 0) with v — is called the 

comoving frame. Transition to another inertial frame are fulfilled by a Lorentz transformation L^n^ = n M . A 
complete set of inertial frames {n^} is obtained by all Lorentz transformations of a comoving frame. 

2 Recall that this specific reference frame was chosen by Wigner jjjj, in order to construct irreducible repre- 
sentations of the Poincare group supposing the existence of a vacuum as a state with minimal energy (see in 
detail piUTj ). 
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These functional were called the "dressed fields", or the "radiation variables" [T]. 

In terms of variables (J0J) - (fTTJ) the Gauss law (pj) takes the diagonal form (i.e. it loses the 
linear term d^A}.) 

AA^ d) (A) = -j'J rad) = -e^ rad ) 7o ^ rad ) (13) 
and the spatial components become transversal identically 

d k A^ d) (A)=0. (14) 

After the substitution of a manifest resolution of the Gauss constraint © into the initial 
action JH) 3 this action can be expressed in terms of the gauge-invariant radiation variables (0) 



'S'QED I *Sqed =q = S'qed - ( 15 ) 



d A x 



l(^4 rad) ) 2 + ^[ip - + 4 rad) j| rad) - ^ rad) ^j'S rad) 



One can see that Dirac in managed to express his Hamiltonian scheme in terms of radiation 
variables © - (fTTj) as functionals from the initial fields, and these functionals are invariant 
(TH^ with respect to gauge transformations of the initial fields JBJ). In other words, one can say 
that the Dirac radiation variables are manifestly gauge-invariant. Therefore, the fact that the 
current jft^ cannot satisfy the conservation constraint 

^j (rad)M + 0. (16) 

does not contradict the gauge-invariance of the initial action JTJ and, moreover, it allows us to 
give the gauge-invariant description of bound states, where the current is not conserved. 



2.3 Dirac 's frame depended approach versus the frame free one 

Values of the reduced action (fl"5|) onto solutions of the motion equations takes the form 

jf ad) (x)C 00 (a; - y)jt d) +3t d \^(x - y)jt^{y) 



QED I ^Q ED -Q ~~ 2 



+ 



+ / rf 4 x^ (rad) [^-m]V (rad) , (17) 



where 



CooOr-y) = -5{x -y )——^ — (18) 

Aix\x — y\ 

T*(* -y) = / A (hi ~ ^f) e ^L. (19) 

(2tt) 4 J \ q 2 J ql + ie 

The action (fl"T|) allows one to determine the propagator of radiation variables in the momentum 
representation j(x) = (l/(27r) 2 ) J d A qe~ %qvxV j (q) . One can see that the current part of the action 
(fTTj) takes the form 

S (iad) [j] = \j d 4 qf(q)f(q)D^(q), (20) 

3 This substitution, i.e. the calculation of value of the action onto a solution of the Gauss constraint, is called 
the reduction procedure. This reduction allows us to eliminate nonphysical pure gauge degrees of freedom [12] . 
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where 

^ (rad) (?) = <W*^ + 8^ (<% - ^) = C, v + T,„ (21) 

is the propagator as the sum of the instantaneous exchange and transverse field one T MJ ,. 

This propagator contains two singularities: the Coulomb one C forming instantaneous 
atoms (see Appendix A) and the light cone ones of the transverse variables T (considered as 
radiation corrections). This propagator can be identically rewritten as the sum of the Feynman 
propagator F (that does not depend on the frame time-axis n M ) and the longitudinal L one: 

s + (gb^)( g o^) - (^)( gj -^) s {q) + ( } _ (22) 

It is easy to see that if the currents ^(g) satisfy the conservation constraint 

g „~(rad) = g (33) 

in contrast to JTBJ), the action (JIQ)) takes the form 

S (rad) b1 = \J d*q3»(q)F„„( q )f(q) = -~ J d 4 q^-. (24) 
This result corresponds to the Lagrangian 

^ED = ^QED-i(^) 2 , (25) 

where the first term is the initial Lagrangian J2J. The change 

£qed|(s m ,7'm = o) £qed = £qed - 2"(^^ M ) 2 ( 2 6) 

violates gauge-invariance (in contrast to the radiation variables) and it is known as the Feynman 
gauge. A similar change (jZEj) is the basis of the accepted frame free formulation of SM 

Thus, in the Dirac Hamiltonian approach to QED we have perturbation theory featuring 
two singularities in the photon propagators (|2T| (the instantaneous singularity and that at 
the light cone), whereas the perturbation theory in the frame free gauge (|26|) (accepted now 
to formulate SM [H]) involves only one singularity in the propagators (that at the light cone 
The perturbation theory in this gauge in terms of the propagators with the light cone 
singularity can by no means describe instantaneous Coulomb atoms, and it contains only the 
frame free Wick-Cutkosky bound states whose spectrum is not observed in nature [T3l fTlj. 

The problem of different physical results in two different gauges cannot be explained by the 
violation of the gauge-invariance in the Lagrangian (|26j). We can keep gauge-invariance of the 
frame free method choosing the Lorentz variables A^\A) by gauge transformations 

A^(A)=A fl + d tl A^(A), ^ L )(^,A) = e ieA(L) ^V, (A (L) (A) = -^d tl A») . (27) 

The gauge-invariant functional Ajt\A) identically satisfies the Lorentz constraint d^A^^ = 0. 
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In fact, the change from the radiation variable (j2J A^ r£ld \ ?/A rad ) to the frame free Lorentz 
ones 

A^(A^) = +d ll A {Tad \A^), (28) 
V/ rad) (^ (L U (L) ) = e ieA(rad) ( A(L) V (L) , (29) 

where A^- Tad ' (A^) is the Dirac phase given by Eq. (JHJ), keeps all physical results, if there remains 
the Dirac gauge of physical sources in terms of the new variables: 

4ourL = + d k A^ d \A^)} + fjWfeA^HAM)^) + ^Jg-^MHAW)^^ (3Q) 

As a relic of the Dirac fundamental quantization, the Dirac phase factors e * eA(rad) (^ (L) ) j n the 
source Lagrangian (|3U1) "remember" the entire body of information about the frame of reference 
to initial data and the instantaneous interaction. As was predicted by Schwinger [16j, all 
these effects disappear, leaving no trace, if these Dirac factors are removed by means of the 
substitution 

^sources ~ * ^sources = j[i ^ ^ + ^ ^ + ■ (31) 

This substitution (treated as another choice of gauge of physical sources by the elimination of 
the Dirac phase factors) loses the instantaneous interaction and instantaneous bound states. 
We see that there is the strong dependence of physical results on the gauge of physical sources, 
and this dependence does not mean the violation of the gauge-invariance principle, because 
both the radiation variables and the Lorentz ones (|27l) are gauge-invariant functionals of 
the initial fields 4 . 

The region of validity of the theorem equivalence of the frame dependent method and the 
frame free Lorentz gauge formulations (changing the gauge of physical sources) is restricted by 
the elementary particle scattering amplitudes |3j. 

The change of the gauge of physical sources (|3T| is just the point where the frame free FP 
method in the Lorentz gauge ((SH) loses the instantaneous interaction (|THI) forming instantaneous 
bound states (considered in Appendixes A and B). Moreover, the field A treated as a dynamic 
one in the the Lorentz gauge formulation (|26|) gives a negative contribution to energy and it 
leads to tachyon in the spectrum of the Wick-Cutkosky bound states p~3| fTi]. 

The standard frame free approach (BJ to SM extending the FP method out of the region 
of its validity finally lost instantaneous interactions. However, if we keep instantaneous 
interactions, the question arises about a choice of the frame and relativistic invariance. What 
do relativistic covariant instantaneous atoms in QED mean? 

2.4 Relativistic invariance: "frame free" versus "many frames" 

It is interesting to trace the historical evolution of the physical concepts of relativistic invariance 
and gauge invariance. Both Julian Schwinger (who "rejected all Lorentz gauge formulations 
as unsuited to the role of providing the fundamental operator quantization" ^H]) and Faddeev 
(who accepted these formulations in Pj) used the same argument of relativistic invariance. The 
phrase "Finally, the Hamiltonian formulation, in the case of field theory, is not relativistically 

4 Therefore, the assertion that "the dependence of physical result on the choice of a gauge means the violation 
of the gauge-invariance principle" is not correct. The correct full name of the concept "gauge" is gauge of 
physical sources. A change of gauge of physical sources can lead to a change of physical results as we have seen 
above. 
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covarianf was treated in [3| as a defect of the frame dependent Hamiltonian formulation and 
one of the arguments in favor of the Lorentz frame free gauge (without the instantaneous in- 
teraction). However, in accord with the theory of representations of the Poincare group the 
relativistic invariance means that a complete set of states {|<E>j >} n d obtained by all Lorentz 



corn- 



transformations of a state |$ >o in a definite inertial frame of reference n° = (1,0,0,0) 
cides with a complete set of states {|$/ >} n obtained by all Lorentz transformations of this state 
|$ >o in another frame of 'reference = (-7=5=, ^7=5=) [IDI- Therefore, the frame dependence 
of the Hamiltonian formulation was treated in (TSJ as necessity to prove that the algebra of 
commutation relations of gauge-invariant observables coincide with the Poincare algebra. This 
proof was fulfilled by Zumino (THj (and Schwinger fUj in the case of non-Abelian fields even in 
1962) at the level of the algebra of commutation relations of the Hamiltonian and momentum 



operators with the nonlocal commutation relation i[A^ ad \x), Pj Tad \y)] = 5kj — S 3 (x — y) 

for quantum radiation variables taken to be symmetrically ordering in the Hamiltonian. 

We can see that the relativistic invariance as the theory of irreducible representation of 
the Poincare group does not mean the frame free formulation IpTjl . The relativistic invariance 
means that a complete set of all physical states includes the states obtained by all Lorentz 
transformations of the comoving frame, i.e. all Lorentz transformations of the instantaneous 
interactions obtained in the comoving frame. Thus, there are two following possibilities to 
construct the weak S-matrix amplitude transitions between physical states, i) the frame free 
S-matrix, and ii) the many-frame S-matrix. 

2.5 Many-frame relativistic S-matrix with instantaneous interactions 

Recall that S-matrix in QED is constructed on the basis of the Schrodinger equation 

# (rad) |$(x°) >= i-^\$(x°) >, (32) 

where 

# (rad) = / d 3 xH^ d) [n]=J d 3 x[no(n)+n- ms : mt Xn)+n TCt . mt .(n)}, (33) 

V Vo 

H {vad \n) = [P/ rad) ,9 4 (rad) + p| ad) <9 ^ (rad) - £ (rad) ] (34) 

are the Hamiltonian and its density as a sum of free field part, instantaneous interaction part, 
and retarded interaction one; $ is a physical states. All these notions (space-time, volume, 
energy, class of functions and etc.) make sense only in a definite frame of reference to initial 

data distinguished by a unit time-axis = ^ v /j~g2 ' > where v is the initial velocity of 

an inertial frame, v = is called a comoving frame (cf) (see the footnote P). The resolution of 
the equation takes the form of the operator of evolution 

|$(a; ) >= Texp ^~ i J o ^°# (md) j |$/(4 n) ) >, (35) 

here ^/(x^) > is an initial state defined in the asymptotic region, where the Hamiltonian of 
retarded interaction in the sum (|33|) can be neglected 5 , so that one can determine the energy 

5 We can neglect the retarded interaction provided that |^( out ) — x (in)l ^Tmin> ^ ^7 mm- This con- 

dition means that all stationary solutions with zero energy Ei^ min — > cannot be considered as perturbational. 



dhdi 
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spectrum Ej from the equation 



d 3 x[H (n) +^ ins , nt .H]|$ / (x° in) ) >= S z |$ r (4 nl ) > 



'(in)/ 



(36) 



In accord with the spectrality principle a complete set of eigenvalues Ej (with a finite energy 
density Ei/Vq < oo) of a manifold of all asymptotic states with the quantum numbers I includes 
a physical vacuum as a state with a minimal energy Ej min 

Notice that in the comoving frame n^ f = (1,0,0,0) the set of one-particle and two-particle 
equations 6 (|36|) coincide with the Schwinger - Dyson equation and the Salpeter one, respectively 
(see Appendix A). 

The next step is the construction of the S-matrix elements as weak probability amplitude 
transitions between the all asymptotic 7 states including bound states. A complete set of these 
states includes any bound state moving with the momentum P M obtained the Lorentz trans- 
formation of the time-axis 



n 



cf 



1,0,0,0) 



v_ T>* 



Therefore, practically, in QED, one uses the many-frame S-matrix 



(37) 



S 



(rad, practical) 

<^out )^in 



A x (out), 



exp 



d A xH { ™ d) (h m ) 



l{ X %a)) 



(38) 



where the time-axis n?f = (1,0,0,0) [T 



tional to the operator of the total momentum 



is replaced by the initial datum operator n^ D propor- 



n™\<S> J >= 



p j 

Mr 



(39) 



of any physical state considered as a irreducible representation of the Poincare group [20} l2T]. 
where this operator n^ f becomes c-number (see Appendix A). 

Thus, instead of the acceptable frame free formulation of the Standard Model without 
instantaneous interactions we suppose to use the Dirac-like radiation variables with instanta- 
neous interactions and the many-frame S-matrix (|3%J) where the comoving time-axis is replaced 
by the initial datum operator (|39|) 9 . 

6 Similar equations are well - known from the nonrelativistic many - body theory (Landau's theory of fermi 
liquids ^2], Random Phase Approximation [Ts] l and play an essential role in the description of elementary 
excitation in atomic nuclei JHj- Their relativistic analogies describe in QCD the Goldstone pion and the 
constituent masses of the light quarks |2()ll21ll^] (see Appendix B). 

7 The hypothesis of asymptotical states of S-matrix in QFT means that physical states are created together 
with their comoving reference frames where quantum numbers of these states are measured. 

8 It is clear that at the point of the existence of the bound state with the definite total momentum Va^ any 
instantaneous interaction with the time-axis parallel to this momentum is much greater than any retarded 
interaction in the sum i-VM \2','<\ . 

9 In particular, this definition of S-matrix IK-S8|1 justifies the application of the axial gauge for the description 
of the deep-inelastic scattering amplitudes, where this gauge imitates the transition to the large momentum 
frame EH- 
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3 Instantaneous weak interactions 



3.1 Standard Model 

We consider the accepted SM without any additional terms violating gauge invariance 



.212 



<P[fv(e R e L + e L eij) + (e R e L + e L e fl )] + 



+^ 2 [9 2 W+W^ + (g 2 + g' 2 )Zl] 



(40) 



£/ = v e %{jd)v e + 



ig sin Qy/A^e + 

sin 2 fl 



zcostV 2cosfcV 



7 w 



(e ltl e)Z^ + ^=(P e W+re L + e L W-r^) (41) 

COSt^y y/2 



Cy = -\(d,A u - d v A,f - l{d„Z v - d u Z,f - \[D,W+ - D V W+] 2 (42) 
- ie{d,A u - d u A,)W + »W-» - g 2 cos 2 9 W [Z 2 (W^W^) - (Wj Z»)(W; Z v )\ + 
+ igt^Q^Z&Z^W^W-^igccnw c.c], 

where V?- = (A^, Z^, W+, PV") are the standard set of vector fields, and with the Higgs effect 
of spontaneous symmetry breaking 



2 m v 

V ' ^ = 2A ' 



(43) 



where rj is the Higgs field with the mass m v . 

First of all, we can see that this action (Jl2*|) is bilinear with respect to the time components 
of the vector fields V K = (A , Z , W+, Wq) in the "comoving frame" rajf = (1, 0, 0, 0) 



d 4 x 



ZV K f KI V I 4- T/ A 4- 



-K T K 



(44) 



where is the matrix of differential operators. Therefore, the Dirac approach to SM can 
be realized. This means that the problems of the reduction and diagonalization of the set of 
the Gauss laws are solvable, and the Poincare algebra of gauge-invariant observables can be 
proved. 

In any case, SM in the lowest order of perturbation theory is reduced to the sum of the 
Abelian massive vector fields, where Dirac's radiation variables was considered in |2*5] . 



3.2 The reduction of the Abelian massive vector theory 

The classical action of massive QED is 

1 „ 1 



W 



d xC(x) 



d x 



(45) 
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with the Lagrangian C(x) and the currents J M = e\l/7 M \l/. The Euler-Lagrange equation for Vq 
is a constraint 

= o <_> (A - M 2 )V = -$V$ + Jo , (46) 

OVq 

It has the solution 

Vo[V, M = X^JP { ~ 9iVi + Jo) • (4?) 
Inserting this into W we obtain the reduced W red 

W rcd [V, *] = J d 4 x£ red = y rf 4 x (jCl d + £* d ) , (48) 



where 



£L = - y^ + ^(A-M 2 )^, , 



= ^o^^Jo + Jo(^^^+^ + *(z^-m)*, (49) 



are the reduced Lagrangians with the operator 



2 



did,; rp M 
where we have used the longitudinal and transverse projection operators 

In contrast to the massless case, R^ is not a projection operator, R 2 ^ R, but 

RuRji = Rji H ■ — ;r (52) 

« " (A-M 2 ) 2 V ; 

jTi-gT A ~ M % H (53) 

In the massless limit however the reduction operator becomes the transverse projection operator 
for photons and ceases to be invertible. 

The second term in Cf ed ijlHJ) can be eliminated (in order to diagonalize the Gauss law) by 
introducing the new radiation (R) variables 

v >" = (** - a^f) v > = ^ (54) 



and Rij is invertible 



* H = exp (ie A _ M2 ^ ) * (55) 

which generalizes the Dirac radiation variables in QED (fTt)|) and (fTTJ) . Since (|55|) is only a 
phase transformation, the corresponding current J M stays the same and Cf ed becomes 

A* d = + jj - A^y ^ + * - m) * R - (56) 
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Using the radiation variables V k R = RkjVj with Rij given by (jHOjl the reduced Lagrangian can 
be written 

£ red = \ (vfR^V* + V*(A - M*)R^V R ) + * R (ip - m)* R + W* + I J ^-L^J . (57) 

Thus, the frame-fixing = (A , A k ) 7 the treatment of A as a classical field, and the Dirac-like 
diagonalization of the Gauss law (jlHI) (i.e. eliminating the linear field contribution diVi) lead 
to radiation variables ijHljl and ijHHJ) that cannot be associated with any gauge constraint, in 
the case of the massive fields [25] . 



3.3 Dirac's frame dependent approach versus the frame free one 

The Fourier transform of the reduction operator 

£ef (A) (?>f (A) (?) = R*{Q) (fo + f|) R l0 {k) = Rij (q) = S kl - -^f- . (58) 

leads to the free propagator 

Dfj(x-y) = (0|T^(f^ )^ B (y,2/o)|0) 

~ l J (2tt) 4 g 2 - M 2 + *e ~ <f + M 2 J ' (59) 

Together with the instantaneous interaction (see the last term in the reduced Lagrangian (foTjl ) 
this leads to the following current-current interaction 



J" J- = Jbjrns. * + J «4 («« - ?fb> ) ifhp ■ (60) 

It is the generalization of the complete radiation photon propagator (|2*T]) in QED. The propa- 
gator (J60|) contains the instantaneous part as a analogue of the Coulomb instantaneous one in 
the propagator (|2~Tf . 

As it was shown in [25], the Lorentz transformations of classical radiation variables coincide 
with the quantum ones and they both (quantum and classical) correspond to the transition 
to another Lorentz reference frame distinguished by another time-axis, where the relativistic 
covariant propagator takes the form 

D > W = W?W\ " ^ " 5^ J (61) 

with = — n^qn), 5^ u = 5^ v — n^n v . Therefore, we shall use the radiation propagator 
(JET}, where n M is treated as the initial datum operator instead of the conventional frame 
free massive vector propagator [E] 

J^(q)r = -j^ J * - W) r (62) 

In contrast to this conventional frame free massive vector propagator the radiation propagator 
(jfinj) is regular in the limit M — > and is well behaved for large momenta. 
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For a better comparison with the conventional covariant propagator (Ifi2l we rewrite the 
propagator (f60|) in the alternative form 

Hence we see that for massive QED, where the vector field is coupled to a conserved current 
= 0), we find that the effective current-current interactions mediated by the propagator 
of the radiation fields (|6T)1) and by the conventional frame free propagator (|fi2l) coincide 

j»D* v r = j^r . (64) 

If the current is not conserved Jo^o 7^ JkQk, the radiation field variables with the propagator 
(fHUf (or ((EH)) is inequivalent to the frame free variables with the propagator (|H2~| because 
different variables corresponds to different gauges of physical sources. 

In the next Section, we show the dependence of physical results on gauge of physical sources. 

4 Observational tests 

4.1 Chiral Bosonization of EW Interaction 

In order to illustrate the inequivalence of Dirac's method (|5T|) and the frame free one [6j, we 
can consider the K ^ ir transition amplitude and the kaon decays, where the dominance of 
the instantaneous interactions was recently predicted Recall that the observation of kaon 
weak decays has been crucial for the modern theory of particle physics fZ7\ . 

It was accepted to describe weak decays in the framework of electroweak (EW) theory at 
the quark QCD level including current vector boson weak interactions [2H1 EH] 

£(J) = -(J»K + J+W-) = — 6 U;W++4W-), (65) 

2 V 2 sin t>w 

where Jt = d'j^l — 75)^; d! = d cos 9c + s sin 9c, 9c is a Cabbibo angle sin^c = 0.223. 

However, a consistent theory at large distances of QCD is not yet constructed up to now. 
Therefore, now the most effective method of analysis of kaon decay physics [30, EH E2 ESI 
is the chiral perturbation theory jSH EE] (the list of arguments in favor of this perturbation 
theory is given in Appendix B). 

The quark content of 7r + and K + mesons n + = (d, u), K + — (s,u),K = (s, d) leads to the 
effective chiral hadron currents Jr in the Lagrangian flHHI) 

4 = [4±*4] cos 9 C + [4±*4] sin 9 C , (66) 
where using the Gell-Mann matrixes \ k one can define the meson current as [3H 

i E A "4 = iXk ( V " - Kt = F*e*%e-*, (67) 



8 

k=l 



1T-V2 

K-V2 



7T 



y/2 K+V2 \ 



(68) 



12 



(b) 



7T°(k + l) 



— » #vVW# — 

K+(k) W + {k) 7i+{k) 



K+(k) 



W + (-l) 



7T 



-(k) 



Figure 1: Axial (a) and vector (b) current contribution into K + — > tt + transition 
In the first orders in mesons one can write 

V- = V2 [sm9 c {K'dy - n°d^K-) + cos6 c (t^tt - n°d^-) ] + ... (69) 

and 

A- = V2 F n (d^K- sm9 c + ^tT cos6 c ) + ...; (70) 
here F n ~ 92 MeV. The right form of chiral Lagrangian of the electromagnetic interaction of 
mesons can be constructed by the covariant derivative d^x^ D^x^ = ipn ± ieA^x^, where 
X ± = K ± ,n ± . 

We suppose also that the quark content of the mesons determines hadronization of QCD 
[2fl| l2T] conserving its chiral and gauge symmetries. 

4.2 K — > 7T Transition Amplitude and the Rule AT = - 

Let us consider the K + — ► tc + transition amplitude 

<n + \-i J dx i dy i J»{x)D™{x-y)J u {y)\K + >= i(2n) 4 5 4 (k - p)G EW i:(k 2 ), (71) 
in the first order of the EW perturbation theory in the Fermi coupling constant 



sin 9 C cos 9 C e _ . G 



GEW = ~~7^stf0 w " ™0c«*0c^ (72) 
comparing two different W-boson field propagators, the radiation (R) propagator (jHT| 



Dl(q) _ w = _ ^ _ j (73) 



where q^ = q^ — n^qn), 5^ v = 5^ u — n^n v in the frame = k^j v/fc^, and the frame free (F) 
propagator (|62|) 



These propagators give the expressions corresponding to the diagrams in Fig[T] 

E(^) E«(P) = 2F X D±(k\n)K + 2. / A ^f^ , (75) 



M ^ V 1 7 " 7 (2tt) 4 (Jfe + qf - ml + ze 

= 2*^(1^ + « / ^ 2k ^ D ^~^ + ^ = (76) 



' " ' 7 (2tt) 4 (fc + g) 2 -m2+ze 

(77) 



2fc 2 ^ + f d ^ (2k lM + q,)(2k u + q u 



"M^ J (27r)*(M^-q2-ie)[(k + qy-ml + ie} 

13 



M 2 W 



Both versions R and F coincide in the case of the axial contribution corresponding to the first 
diagram in Fig. [Hand both they reduce to the instantaneous interaction contribution, because 



k^D F (k) = hn?D%(k\n) 



M 2 W - 



However, in the case of the vector contribution corresponding to the second diagram in Fig. [T] 
the radiation version differs from the frame free one (fTTj) 10 . 

In contrast to the frame free version (fTTf) . two radiation variable diagrams shown in Fig. [T] 
in the frame ~ = (ko, 0, 0, 0) are reduced to the instantaneous interaction contribution 



i(2n) 4 5 4 (k - p)G EW E R (k 2 ) =< 7T + | -i I dx*J (x) 



Ux)\K + > 



(78) 



with the normal ordering of the pion fields, so that 

Fl 1 



E R (k 2 ' 



2k 



Ml 



d 3 l 



(79) 



here E n (l) = y m\ + P is the energy of 7r-meson. The reduction of the radiation variable 
loop diagrams to the instantaneous interaction ones demonstrates that the radiation variables 
allow to separate the low-energy region from the high-energy one. This separation justifies 
the application of the low-energy chiral perturbation theory {36 j . where the summation of the 
chiral series can be considered here as the meson form factors [SOIESIES- Finally, the radiation 
variables (|THJ) and (f78|) give the result 



E R (k 2 



2k 



M 2 W 



+ 



d 3 l fKww(—l )f; 



irnW 



-I 



7^ 



(2tt) 3 7 2E w (f) 



M^ + P 



2k 



Ivl w 



9s; 



(80) 



here fx^wi^ 2 ) an d f^wi~ ^ 2 ) are ^ ne wea k vector form factors introduced in the chiral per- 
turbation theory |36| as the natural regularization of the low-energy meson instantaneous in- 
teractions by the next orders, and 



28 = 1 + 



M 2 W 



(2tt) 2 F 2 J e^{1) 



d\i\ i 2 fLwH 2 )flwH 



7^ 



M^ + P 



(81) 



is the parameter of the enhancement of the probability of the axial K + — ► n + transition. The 
relation JEH), where g$ = 5.1 {SIlEIj, can be considered as the low-energy sum rule for meson 
form factors [38J in the space-like region. 

This result shows that the instantaneous vector interaction can increase the axial interaction 
K + — > 7r + transition in g 8 times and give a new term describing the K° — > n° transition 
proportional to g% — 1. 

Using the covariant perturbation theory [39J developed as the series J„(7 © £) = J_u(Q + 
F 2 d^ k + 7* fijkJ_i(£) + 0{pf 2 ) with respect to quantum fields 7 added to £ as the product 

10 The Faddeev equivalence theorem [2j is not valid, because the vector current J M = Kd^n — ird^K becomes 
the vertex — Kd^D^ — D^d^K, where one of fields is replaced by its propagator DD^ = 5(x), and d^T 11 ^ 0. 
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)7 *(g) ( & ) <Y(q) 

■* — #^v#— — i>— — — > — m—* — */vw » 

W+(jfe) 7T+(A;) tt+Qo) AT+(fc) tf+(p) PF+(p) 7r+(p) 




Figure 2: Axial current contribution 



e *7 e «£ = e *(7©C) [3gJ 5 one can see that the normal ordering N(x — y) =< 0\Y(x)Y' (y) |0 >- 
<5 M ' J d 3 le d '( x ~v> /(2E n (l)) in the product of the currents can lead to an effective Lagrangian 



<3 _N(z)e- Mw ^ 



4n\z\ 



JL{t{x))£(Z{z + x))[(f m + if ij2 ){f m - if if5 ) + h.c] < 0| 7 WQ/)|0 > • 



In the limit Mw — > oo, one can neglect the dependence of the current on z, and we get the 
effective Lagrangians [471] 



Gf 

£(At=±) = ~^98 cos 9 C sin 9 C 



1 



(4 + *4) (4 - *4) - (4 + 7f 4) 05 - *4) + h - c - > ( 82 ) 



CtAT-3\ = — 7= cos sin #c 



1 



(4 + -^4)(4-i4) + ^. c . 



(83) 



This Lagrangian at the level of the tree diagrams describes the nonleptonic decays in satisfactory 
agreement with experimental data within the accuracy 20 -j- 30% [34") 14*0]. 



4.3 The K+ -> tt+ + Z + / amplitude 

The result of calculation of the axial contributions to the amplitude of the process K + — > 
7r + + I + / within the framework of chiral Lagrangian 1)650 . ([Hfijl including phenomenological 
meson form factors shown in Fig. [2] as a fat dot takes the form 



T; 



AA 

(K+^TT + l+l-) 



2eG EW L v D^ ad \q)(k,+ PlM ) t AA {q\k 



2 2\ 
,P ), 



where Gew is the constant ((72*|) . L M = Ij^l is leptonic current 



(84) 



(85) 



and (? 2 ) are meson form factors. In this expression the propagator Dj} u (p\n) given by Eq. 



Tj keeps only the covariant instantaneous part 
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7T 



'(* + /) 



>7*(?) 



+ — ^w^- 



K + (k) W + {~1) Ti + {k) TT+Qo) 



(6) 



5 7*(g) 7T°(p + /) 



fT + (£;) tf+(p) W + (-Z) 7r+(p) 




W+(-Z) 7r+(p) 



(d) 7 *( ? )<^ 7T°(p+Z) 




W+(-Z) 7r+(p) 



Figure 3: Vector current contribution 



On the mass-shell the sum l|%H|) takes the form 



A'- 



% 2 



eVi 



mt 



M\ — ml 



(86) 



As can be seen in (HOI EH- the amplitude for K + — > 7r + + Z + Z vanishes at the tree level, 
where form factors are equal to unit t(q 2 )\fV = fA =l = 0. It is well known that these form factors 



in the limit q 2 — > take the form f^'iciQ 



< > 



TT,K 1 



where < r 2 > 



(A,V) 



are the axial 



and vector mean square radii of mesons, respectively, in agreement with the chiral perturbation 
theory [HI El EE] • 

The result of calculation of vector contributions to the amplitude of the process K + — »■ 
7r + + Z + Z within the framework of chiral Lagrangian (|65j). ()12ft including phenomenological 
meson form factors shown in Fig. El as a fat dot takes the form 

„2\ 



1 (K + ^TT + l+l-) 



2eG EW L v D^ ad \q)(k,+ Pfl ) t vv (q 2 ,k 2 ,p 2 
where we use the propagator (|6T1) in the initial datum frame (|3H 



(87) 



t vv (q 2 ,k 2 y 

fW)p 2 

M 2 K -p 2 - ie 



2(2n) 2 J En (l) 



mt 



k 2 — ie 



+ 



here /^^(Z 2 ), /avha^Z 2 ) are the four particle interaction form factors which should coincide 



with the three particle interaction form factors f^yy-y 
with the gauge invariance of the strong interactions. 



■kttWi J KnW-y 



fif-ur in accordance 



Using the result (JHH) one can write on the mass shell k 2 = Mj^p 2 



m 



KnW 
2 . 



t vv (q 2 ,M 2 K ,ml) = (gs-mq 2 )., 
t AA (q 2 ,M 2 K ,ml)+t VV (q 2 ,M 2 K ,mi) = g 8 t(q 2 ), 



(89) 
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where gs,t(q 2 ) are given by Eqs. (j£Tl) and (JHEJ), respectively. 

Finally, the result of calculation of axial and vector contributions to the amplitude of 
the process K + — > tt + + 1 + 1 within the framework of chiral Lagrangian (jl|,(pD including 
phenomenological meson form factors takes the form 

T^l n+l+l _ } = 2eg 8 G EW L iy D^ ad \q)(k, + p ll ) t(q 2 ). (90) 
This amplitude leads to the total decay rate for the transition K + — > 7r + e + e~ 



2 



where 



- da 2 

T iK+ ^ n+e+e ^((M K -m T ) 2 ) =T e+e - I ^p( g 2 )F(g 2 ), (91) 



Ml 



- (s lCl c 3 ) 2 g 2 G 2 F a 2 M% n7yir 22 rpV 



F{q 2 



and 



{A-K)H{q 2 ) 



V(2\ 



eV f„2\ 



mi 



'M 2 K -m\ 



(93) 



4m 



2\ 1/2 



x 1 + 



2mi 



q L 



A 3 / 2 (l,g 2 /M 2 K ,ml/ 'M 



here \(a,b, c) = a 2 + 6 2 + c 2 — 2(a6 + 6c + ca), S1C1C3 is the product of Cabibbo-Kobayashi- 
Maskawa matrix elements V u dV us . 

The mechanism of enhancement |AT| = | considered in this paper can be generalized to a 
description of other processes including K + — ► 7r + z/P by replacing the 7-propagator by the Z- 
boson one, so that we have the relations 



9 Z 



(47Tig 2 J p(g2 )r 



Mi 



<ir( A -+_ 7r + e + e -)(g 2 ) 



M 2 Z 



(4ttF w ) 



/i(<7 2 )+/.V 



M| rfT(A:+->^+ y i/)(g 2 ) _ 
p(g 2 )r, p rfg 2 

■ £V) + tw) - /:(^ 2 )] 



M 



-1 2 



(94) 



Thus, exploiting the instantaneous weak interaction mechanism of enhancement in the K — > 7r 
transition probability and QCD symmetry we derive the sum rule of EW vector meson form 
factors given by Eqs. (f8Tf . (JHEJ), (|9Tf and their relation to the differential K — > ne + e~ decay 
rate (JHU) - 
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4.4 The form factor probe of the differential K — > ire + e decay rate 

The estimation of the meson loop contribution was made in [ST] where a function 4> 2 {q 2 ) was 
used instead of the form factor rate F(q 2 ) ("9""| 

F(q 2 ) => 2 (g 2 ) 

It was shown that the values g$ = 5.1, 0(0) = 1 gave the total decay rate T( K +^ w + e + e -) = 
1.91 x 10~ 23 GeV in the satisfactory agreement with the experimental data T^+^Tr+e+e-) — 
1.44 ± 0.27 x 10~ 23 GeV (S3 SI]- However, the main contribution goes from the baryon loops 
"SHIES ESI- Therefore, we discuss here the value of the differential K — > ne + e~ decay rate ("9""f 
in the chiral perturbation theory [SU EB] where both the pion loop contribution and baryon 
ones lead to the meson form factors and resonances [SDl EE] in the Pade-type approximation 

1 



f v Aq 2 ) = 1 + M p - 2 g 2 + aoU^q 2 ) + ...- - ^-f - ^ (95) 

1 



where 



/.V) * f^q 2 ) = l + M- 2 q 2 + ...^-— — — (96) 

^ -'"a q 

ml 4 2 

°" = (4^3 = 103' (97) 

f\ \ 1/2 / t 1 / 2 \ 

U n (t) = (1 - t) {■= - lj arctan ( — — j- f - j - 1; < t < (2m,) 2 (98) 

^ t-ir i\ 1/2 r , W 2 + a- iW 2 ) , , 2 
n,(t) = — (^i - 1 J (^-iog t - 1/2 _| t -_/ )1/2 j + i; (K) 2 <t (99) 

is the pion loop contribution [SUES], and M p = 771 MeV and M a = 980 MeV are the values 
of resonance masses [Sl| that can be determined by the baryon loops [301 E2 ESI EE] ■ 
The constant approximation of the type of [31 J 



F{q 2 = 0) = F(0) = (4nF n f [f ^ q2) f J {q2)]2 \ g , =0 = {A<kF^[M; 2 - M~ 2 ] 2 = 0.74 (100) 

corresponds to the value of TfK+-^+e+e-) = 1-41 x 10 -23 GeV in satisfactory agreement with 
the experimental data T {K +^+ e+e -) = 1.44 ± 0.27 x 10~ 23 GeV |37l PTj . 

The form factors ([9*51) and (T9~""| below the two particle threshold q 2 < 4m 2 determines the 
differential rate F(q 2 ) as a function of q 2 

(101) 

At the region above the two particle threshold 4m 2 < q 2 < (M K — m^) 2 there is a jump of the 
differential rate F(q 2 ) at the level of 5 4- 10 % 

p/ 2\ i a p J M ; 2 - M a 2 f + gfrV ~ ^ 2 )/(4g 6 ) 

F(g } = ^ [1-M-V] 2 [1-M-V] 2 (1 ° 2) 

These results ()100|) . (jlOlj) . and ()102|) are arguments that the detailed investigation of the 
differential K — > Tre + e~ decay rate in the NA48/2 CERN experiment can give us information 
about the meson form factors. To make much more realistic analysis, we have to use the 
unitary and analytic model for meson form factors [38J describing very well all experimental 
data. 
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5 Conclusion 



The instantaneous interactions are inevitable consequence of the Dirac approach to gauge 
theories. Recall that Dirac eliminated all fields with zero momenta (and possible negative 
contributions into the energy of the system) by resolution of the Gauss constraint and its diag- 
onalization, in order to have a physical vacuum as a state with minimal energy in a comoving 
inertial frame. 

This elimination unambiguously leads to the radiation variables associated in QED with 
the Coulomb gauge constraint. However, in SM, the radiation variables of the massive vector 
bosons do not correspond to any whatsoever gauge constraint. The transition from the radia- 
tion variables with the instantaneous interactions to a frame free formulation [6j without the 
instantaneous interactions is fulfilled by gauge of physical sources. The change of the gauge 
of physical sources is just the point where the accepted frame free methods [6j lose all instan- 
taneous interactions together with their physical effects, including the instantaneous bound 
states. The strong dependence of physical results on the gauge of physical sources does not 
mean the violation of the gauge-invariance principle. This means that frame free method has 
a region of validity restricted by the initial datum free processes of the type of the elementary 
particles scattering ones [3], as it follows from the physical meaning of the concept of the frame 
revealed by its full name frame of reference to initial data. 

In order to keep relativistic invariance of S-matrix elements between the instantaneous 
bound states, the instantaneous interactions are determined in an arbitrary frame, where its 
time-axis is treated as the operator of initial data acting in the complete set of the all asymp- 
totic physical states including the bound states considered as irreducible representations of the 
Poincare group [9J. 

In the paper, we compare this Dirac operator approach to SM with the acceptable frame 
free method [6J using as example the K + — ► 7r + transition and the semileptonic and nonleptonic 
kaon decay probability amplitudes. This comparison shows us that the Dirac approach to SM 
separates the low-energy contributions from the high-energy ones and expresses these ampli- 
tudes in terms of the hadron electroweak form factors. Therefore, the obtained amplitudes 
allow to extract information about form factors of ir and K mesons from the K + — > 7r + transi- 
tion in the K + — > n + e + e~ (n + yT , vv) and K + n~ — > e + e~(n + n~ , vv) processes. The parameters 
of the differential K — > 7re + e~ decay rate in the NA48/2 CERN experiment were estimated in 
the chiral perturbation theory. 

These results show that the Dirac formulation of Standard Model can reveal new physical 
effects in comparison with the frame free formulation |HJ used now for describing observational 
data. Therefore, the problem of the Dirac formulation of SM becomes topical in the light of 
future precision experiments. 
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Appendix A: S-matrix for bound states in electrodynamics 

A.l. Spectral problem 

Let us keep in the QED action only the instantaneous interaction (fT%|) neglecting radiation correc- 
tions 

Sqed^H / d 4 x[^(x)(ip-m)^x) + ^ j d*y&(yMx))K(x,y\n)&(x)$(y))}; (A.l) 

here 

JC(x, y\n)= y4 d y(z J -)<5(z • n cf )fi cf , (fi = = 7 • n) (A.2) 

is the instantaneous kernel depending on the relative coordinates z = x — y and Viz- 1 ) = — - — : — r- 

47r|2:- L | 

depends on the transverse components zj^ = z^ — n c ^(z • n cf ) of the relative coordinate with respect to 
the time axis n°* = (1,0,0,0). 

It seems that a most straightforward way for constructing S-matrix of bound states in gauge theories 
is the redefinition of action (|A.1|) in terms of bilocal fields by means of the Legendre transformation 
[42j of the current - current interaction with a kernel IC(x,y\n) 

~ J d 4 xd 4 y^(y)^(x))IC(x,y\n)^(x)^(y)) = (A.3) 
- / d 4 xd 4 yM(x,y)IC- 1 (x,y\n)M{x,y) + [ d 4 xd 4 y{i;(x)^(y))M(x, y). 



2 _ 

and the Gauss functional integral over the bilocal fields in the generating functional 

^(^-(rad^y DMD ^ D $ exM^, (-G 1 + M )) -| (M , XT 1 M ) + i(^ rad ) )+(?? (rad) v| ; (A.4) 
here we used the short - hand notation 

d 4 xd 4 y^{y)i){x){ip-m°)5 {4) {x-y) = (^,-Gq 1 ), (A.5) 
d 4 xd 4 y(il){x)i){y))M{x,y) = (iH>,M). (A.6) 
The integration over the fermi fields tp, rp in l)A.4|) leads to the generating functional 

Z FP [^ iad \fj^] = J DMex V {iW eff [M] + (fj^v^A-Go 1 +M)" 1 )} (A.7) 
with the effective action 142 



W eff [M] = ~(M,lC- l M) + iTrlog(-G 1 +M). (A.8) 
The first step to the quantization of the effective action is the determination of its minimum 



SWeffjM) = _ K -! 

5M ' G^-M 



K ~ M + — 77lM=E = (A.9) 

o 



known as the Schwinger - Dyson equation of the fermion Green function 
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we denote the corresponding classical solution for the bilocal field by = / d i qY^(q ± )e iq ' z ^ S(z ■ n). 
The next step is the expansion of the effective action around the point of minimum M. = S + M., 



Weffp + M) = W^f f + W in f, 

t \7\ — w 

eff 



^{mx- 1 m\ + 1 -(gzM) 2 



oo oo 1 

W mt = Y J W^ = iJ2~(G*M) n , (Gx = (G 1 -Z)- 1 ), 

L — ' L — 4 n 



(A.ll) 
(A.12) 

(A.13) 



n=3 



n=3 



where the small fluctuations of the bilocal field M.(x, y) = Ai(z = x — y\X = (x + y)/2) in the effective 
action (jA.ll|) can be presented as the series of a complete set of the solutions T of the classical equation 
(known as the Bethe - Salpeter one) 



5 2 W eff (Z + M) t 
5M 2 



*M=0 



• r = o 



(A.14) 



with a set of quantum numbers (A) including masses M A = y/'Pfi, and energies u>a = y'P 2 + M\. In 
the momentum representation this series takes the form 



M{z\X)=Y J f 



d?V A^V iP-A,. , i, n , +( .. n , , .-,T,v, 

( (2vr)3/2^2^ (2vr)4 



-{^■ x T A { q ^\V)a + A (V) +e- iV - x r A (q ± \ -V)a A (V)}, (A.15) 



where Ai(z\X) ~ 5(z ■ n) and 

^ A ' ^ 'A! 

the bound state creation and annihilation operators that 

obey the commutation relations 



a AI (V'),a+(V) 



8 A 'aS 3 (V'-V) 







In the momentum representation T(k , V), Eq. (|A.14|1 takes the form 



T(k ± ,V)=i 



d 4 c 



(27T) 



V{k L -q L )ii G Sa (q + V/2)T(q ± ,V)G j:b (q-V/2) ft, 



(A.16) 



(A.17) 



where Gs a , Gsb are the fermion propagators of the two particles (a) and (b). Solutions of equation 
1|A,14|) satisfy the normalization condition ^3] 



d f d 4 q 
i-rzr I ,„ r-, tr 



The corresponding Green function of the bilocal field takes form 



2uja 



(A.18) 



G{q > P |P) = Vl V>o-»H-ie)2»H (V -. H - l e)2. H j ' (A ' 19) 

The Bethe - Salpeter field T{q L ,V) is connected with the wave function by the relation (see [27)| \21\\ 

'V ■ dq 



i 

2^ 



M A 



G Sa (g + P/2)r(g^,P)G s6 (g-P/2) 



(A.20) 



In the nonrelativistic limit in the frame n = = (1,0,0,0) this wave function reduces to the 

Schrodinger one (|A.22|1 

*?®M A ) = (^y^M^ (A.21) 
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satisfying the standard Schrodinger equation 

("27^ + V{ * ) ) - Sch{?) = - Sc/l(f) (A - 22) 

where ^_s c h(z) = / e^ z ^ Sch(o) is a normalizable wave function (f d 3 z\\ ^_sch( z ) l| 2 = 1) in a 

comoving frame = (Ma, 0,0,0); here Ma = (M p + m e — e) is the mass of an atom, M p ,m e are 
masses of proton and electron, /i = M p ■ m e /(M p + m e ) is the reduced mass, and Zk = (x — y)k are 
relative coordinates. 



A. 2. S-matrix elements 



The next step is the construction of the relativistic covariant S-matrix elements considered as 
probability amplitude transitions between the all asymptotic states. These states includes bound 
states given by the spectral equation and all Lorentz transformations of the comoving frame 
distinguished by the time-axis . The Lorentz transformations of the time-axis — > ra( ID ) mean 
the changes of the initial data (ID) of the bound states, i.e. their total momenta Va\i = M^n^ ID \ 
These bound states correspond to the instantaneous interaction with a new time-axis n^ ID ^ 

K{x,y I n( ID )) = fi^V( Zfl - n^z ■ n^)S(z ■ n^)^ 1D \ (A.23) 

One can see that in the space of the bilocal fields 1|A.15|) this unit time-axis n^ ID ^ is proportional to an 
eigenvalue of the total momentum "operator" (i.e. the derivative with respect to total coordinate) 

n^M(z\X) ~ lJL M {z\X) (A.24) 

in agreement with the Markov - Yukawa definition of the relativistic simultaneity in terms of bilocal 
fields M(z\X) = M(x,y) [43J 



-MT^-^i V ,: (A .25) 



Recall that this Markov - Yukawa definition 03] of the relativistic simultaneity in terms of bilocal 
fields has a deeper mathematical meaning as a constraint of irreducible nonlocal representations of the 
Poincare group [H] for an arbitrary bilocal field Ai(x,y) = M(z\X) [2771 144j . 

The relativistic covariant constraint-shell quantization gives both the spectrum of bound states 
and their S-matrix elements. 

Using the decomposition over the bound state quantum numbers (H) 

^H(ab)(Q X \V) = Gxa{q + V/2)T H{ah) {q^\V), (A.26) 

we can write the S-matrix elements for the interaction W^ n ' (| A . 1 3|) between the vacuum and n bound 

states 



< H{P^...,H n V n \iW^ >= -i(2vr) 4 5 4 J>i 

\i=i / j=i 



1 



1/2 

M^(V!,...,V n ) (A.27) 



{ik} 



,-r,".. -"I / 2(^2 + •■• "Pin-l) "Pil 'Pin isp 
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where {ik} denotes permutations over i k ; the amplitude < H{P^ HiVi\iW^\Hi +1 Vi +1 , ...,H n P n > 
can be obtained by the change of the momentum signs Vi+x, ■■■,V n - » —'Pi+i, ■••> —T^n- 

Expressions (|A. . (|A.19|I . 1)A.26|) . I|A.27|) represents Feynman rules for the construction of a 
quantum field theory with the action (|A.11|I in terms of bilocal fields. 

It was shown ^Hj that this quantum field theory with a separable potential of the instantaneous 
interaction leads to the well-known Nambu Jona-Lasinio model HB1 EHj and the phenomenological 
chiral Lagrangians [36j used for the description of the low-energy meson physics. 

In the context of the Dirac approach to gauge theory, to solve the problem of hadronization in 
QCD, one needs to answer the following questions: i) What is the origin of the "separable potential" 
of hadronization in the non-Abelian theory? ii) How to combine the Schrodinger equation for heavy 
quarkonia (that is derived by the residuum of poles of the quark Green functions) with the quark 
confinement? iii) What is the origin of the additional mass of the ninth pseudoscalar meson? 

Appendix B: Instantaneous interactions in QCD 

B.l. QCD action and monopole vacuum 

QCD was proposed [50J as the non-Abelian SU C (3) theory with the action functional 

W = f d^xl—F^+^^ + A^-m]^ , (B.l) 

where ip and $ are the fermion quark fields, = g^A® , 

*g = d A* - Df(A)A b , Ff k = d s A% - d k A] + gf ahc A)A% = e ljk Bf (B.2) 
are non-Abelian gluon electric and magnetic tensions, and 

Df(A) = 5 ab di + gf acb Ai (B.3) 
is the covariant derivative. The action (jB.l|) is invariant with respect to gauge transformations u(t, x) 

A^(xo,x) = u(x ,x) (^A^ + u' 1 (x ,x), ip u = u(x , x)ip . (B.4) 

The class of functions with finite energy density includes the monopole - type functions A$ ~ 0(l/r) 
like the Coulomb potential. 

The asymptotic freedom phenomenon is considered as one of the arguments in favor of that the 
standard perturbation theory in QCD is unstable (SHIES]- One of candidates of the stable vacua in 
the Minkowskian space-time is the stationary monopole - type solution of the classical equation 

A v ~(xo,x) = 0, Al a ~(x ,x) = 6 fc (0,f) = -i^fe iak ^f(r), (B.5) 

where the antisymmetric SU(3) matrices are denoted by := A 2 , \\ := A 5 , \\ := A 7 , and 
r = \x\, the vacuum field contains only one function f(r). 

In the following we discuss the SU(2) case A^ — ► r a , (a = 1,2,3) In this case, the classical 

equation for the function / in Eq. (|B.5p takes the form 

j2 f f/f2_i\ 

Dt h mF b kj ($ i )=Q * -^ + m ^ = Q- (B-6) 

We can see three solutions of this equation fi = 0, f% = +1, fz = — 1 (r / 0). The first solution 
corresponds to the naive unstable perturbation theory. Two nontrivial solutions are the Wu-Yang 
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monopoles [HZ| applied to the construction of physical variables in [SB]- The Wu-Yang monopole is 
a solution of the classical equations everywhere besides the origin of the coordinates r = 0. The 
corresponding magnetic field is Bf($>^) = x a x % /(gr 4 ). Following Wu and Yang j^Zj, we consider the 
whole finite space volume, excluding an e-region around the singular point. To remove a singularity at 
the origin of coordinates and regularize its energy, the Wu-Yang monopole is considered as an infinite 
volume limit (Vq — > oo) of the Bogomol'nyi-Prasad-Sommerfield (BPS) monopole (HHl 



/ 



WY 



1 



/ 



BPS 



1 



(B.7) 



e sinh(r/e) 

with the finite energy J d 3 x[B?(3> k )] 2 = V < B 2 >= A-n:/(g 2 e) = l/(a s e) determining the finite 

volume regularization of the monopole size e = [a s Vo < B 2 >] _1 . The vacuum energy-density of the 
monopole solution < B 2 > is removed by a finite counter-term in the Lagrangian. 



The color particle excitations .4^ rad ^ 

i(rad) 



A (x ,x) = A y Q M, (x ,x), A k (x ,x) = $ fc (ac) + A k ra,d \x , x). 



(B.8) 



are considered as weak deviations in the vacuum background l|B.5|) by analogy with the Dirac radiation 
variables in QED JEJ, d) with the constraint Df rad) (A) = and the Gauss law 



A afc (cD)4 (rad) 



.a(rad) 

~3to ) 



(B.9) 



where A ab (<l>) = [D 2 (&k)] ab is the Laplacian in the vacuum background, and j" is the total color 
current. 

The action can be expressed in terms of the radiation - type variables |2Tj 



(B.10) 



s Q cd|^ =0 - = J d*x [i^Gg" 1 ^ +4 rad) if d) 

+ (V) ( r ad)^(rad iGm) _. 



dx°d 3 xd 3 yj^\x)v ab (x - y) 3 T A \y) + ... + SsocXv^,^) 



■6(rad) . 



where 



s,„c.(»f ( ™ d ,>) <r * d ; 



fB.lll 



is the source action, G^ b and G m are the Green functions of the gluon and quark in the background 



field 



V ab {z) 



z a (z)z b (z) 



^ + E e a a (z)e b a (z)(d.\3\ l -+d + \^ 



a=l,2 



\z\ 2 Att\z\ 

is the potential [22 EC] considered as the Green function 

[A ab (^{z))]V bc (z) = 5 ac 5 3 (z); 



(B.12) 



here z = x — y, e%{z) are a set of orthogonal vectors in color space z a e^(z) = 0, d± are constants, and 
L|_, I- can be found as roots of the "gold section" equation I 2 + I = 1: 



1 + y/5 



-1.618 



L 



-1 + V5 



0.618 . 



(B.13) 



The solution of the set of Schwinger - Dyson and Bethe - Salpeter equations for the rising potential 
was considered in numerous papers, where the effect of spontaneous chiral symmetry breaking was 
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described (see [201 EH EH)- Therefore, this Green function contains the rising potential as the origin 
of the 'hadronization' of quarks and gluons in QCD [61J. 

The Wigner - Markov - Yukawa relativistic generalization of the instantaneous interaction in this 
case can lead to the bilocal theory with the spontaneous chiral symmetry breaking in the low-energy 
limit with the effective chiral Lagrangian one [HUES]- 



B.2. Zero mode of Gauss' constraint 



The BPS - type regularization (|B.7|) of the monopole vacuum is interesting by the existence of a 
nonzero solution of homogeneous equation 



[A ab ($)]Z a = 



of the type of the BPS solution jHS] 

Z(x) 



-ITX 



— /<f PS (0, fo BPS (r) 



1 



(B.14) 



(B.15) 



tanh(r/e) 

The complete solution of the Gauss equation l|B.9p . in this case, is the sum of a solution of the 



(B.16) 



homogeneous equation l)B.14|l and a solution of the nonhomogeneous equation: 

^ (rad) (x) = c(x )Z" - J d 3 yV bc (x-y)$ Tad \y), 

where c(xq) is the zero mode variable. 

The BPS vacuum model gives the possibility to compare the radiation variables in the action (|B. 1 
with the set of the Dirac - type variables removing the field Aq by the gauge transformations 



(D) 

k 

(D) 
I 



A 



(B.17) 



where 



[/(D) 



exp 



dx°c(x°)Z(x) } = exp{N(x")Z{x)}. 



(B.lf 



The Dirac variables l|B.17|l (dressed by the zero mode phase factor) differ from the radiation ones in 



the action QB.lOjl . The differences are Abelian U(l) anomaly and the new gauge of physical sources 



(D)r,(D) 



(B.19) 



5soc.(r/ (D ,r? (D )= / d A x 

as consequences of the gauge transformations (|B.17|) • 

The zero mode of the Gauss constraint l|B.14|) leads to the vacuum electric tension lEBj 

*Svbc = c(x )Df($)Z b = N D?($)Z h (N = d N). (B.20) 

The BPS monopole electric tension is proportional to the vacuum magnetic one. If we choose the 
coefficient of the proportionality as 



Df ($)Z b = Bf ($) 



or 



1 



(2tt) 2 Vq < B 2 > 



47T 



a. 



(B.21) 
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the zero mode c(xo) = N can be identified with a variable given by the winding number functional 



v 



9 C dx° [ d 3 xF"F a » u 



167T 2 / — I t* u 



g | d 3 xDf($)Z b Bf($) J dx°N = N(x° out )-N(x° n ) 



showing that N can interpolate between different classical vacua like an instanton |51| n 

A^out) = "out, ^V«) = ni n . (B.22) 



The equations 1|B.20|1 and 1|B.21)1 lead to the vacuum action of the zero mode dynamics 

W QCD - jrtj - J (|) * . (B.23) 

V 

As we pointed out above, one of consequences of the dressed variables (|B.17|) is the appearance of the 
Abelian anomaly 0H| in the effective meson action obtained by the bosonization of QCD of the type 

Weff = J dx °{\ (% 2 " ™W ) Vb + C^nX (B.24) 

in the pseudoscalar tjq meson isosinglet channel, where C„ = -fy-, Nf = 3 is the number of flavors, 
and F n = 92 MeV is the weak decay coupling. This action is known as the Veneziano one |(i4j . 



The diagonalization of the total Lagrangian supplemented by the action 1|B,23|) leads to an ad- 
ditional mass of the pseudoscalar r/o meson Am, 2 = —J- s ^ — — (2^- This result allows us to 

pz 27T 2 

estimate the value of the vacuum chromomagnetic field in QCD [HE] 

2 _ 2^ 3 F 2 Am r? 2 _ OMGeV 4 

After calculation we can remove infrared regularization Vo — » oo. 

Thus, there are observational arguments in favor of that physical fields are identified with the 
dressed one (|B . 17|) . However, the gauge of dressed sources (jB.19|) differ from the radiation ones (|B . 1 lj) 



by the phase factors ()B.18fl depending on the winding number variable N, The initial position of 
winding number variable A" are degenerated like the U(l) initial data e tX( - 1 ^ = l,xr in \ = 2irk, where k 

is an integer number. The phase factors l|B.18|) are determined within the matrixes v^ n \x) = e nZ \ x ) , 
where n is an integer number which counts how many times a three-dimensional path v^ n \x) turns 
around the 5f7(3)-manifold when the coordinate X{ runs over the space where it is defined. Even, if we 
choose the zero electric monopole energy fixing the zero momentum Pjy = 0, Nr^ = 0, we do not know 
a position of winding number variable due to this topological degeneration of all physical states with 
respect to the phase factor l|B.17|l = e nZ ^ x \n = 0, ±1, ±2, .... Therefore, the creation amplitude 
of the dressed color particles ip( D ^ = ip( Tad W D should be averaged over all degenerated states 

n=+L ^^(in) 1 n =+L 

— V < vacuumlV' (rad) ^ D |g >= 6 , — V e ni ^ = 0, (B.25) 
2L ^ r W |y (2vr) 3 / 2 ^2L ^ r K ' 

n=—L \ / * n= —L 



11 In this case, the monopole interpolation amplitude eKp{iPN[N(x® ut )—N(xf n )]} coincides with the instanton 

g 2 g 7I -2 

one exp{ r-(n ut — ^in)} for a nonphysical value of the momentum Pjy = * — 5-. 

T 9 
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where \q > is the Dirac one-particle state (|B. 17|l with a quantum number q. Averaging over all 
parameters of the topological initial data can lead to a complete destructive interference of all color 
amplitudes [H2] and the zero probability amplitude of creation of color quarks or gluons; whereas 
their radiation propagators have poles like the propagators of electron and photons and the colorless 
amplitudes of the type of the current production < jj^ a j^ a >=< j^ rad ) a j( rad ) a > do not disappear. In 
this case, only colorless ("hadron") states have to form a complete set of physical states. 

Disappearance of color physical states due to the topological degeneration (|B.25|I can be considered 
as color confinement, compatible with the quark - hadron duality as the accepted method of analysis 
of observational data |24| . 
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